LOCAL HOLDER CONTINUITY FOR DOUBLY 
NONLINEAR PARABOLIC EQUATIONS 



TUOMO KUUSI, JUHANA SILJANDER AND JOSE MIGUEL URBANO 

Abstract. We give a proof of the Holder continuity of weak solutions 
of certain degenerate doubly nonlinear parabolic equations in measure 
spaces. We only assume the measure to be a doubling non-trivial Borel 
measure which supports a Poincare inequality. The proof discriminates 
between large scales, for which a Harnack inequality is used, and small 
scales, that require intrinsic scaling methods. 



We consider the regularity issue for nonnegative weak solutions of the 
doubly nonlinear parabolic equation 



This equation is a prototype of a parabolic equation of p-Laplacian type. Its 
solutions can be scaled by nonnegative factors, but in general we cannot add 
a constant to a solution so that the resulting function would be a solution 
to the same equation. 

The purpose of this paper is to obtain a clear and transparent proof for 
the local Holder continuity of nonnegative weak solutions of (jl.ip . Our work 
is a continuation to [17j . where Harnack's inequality for the same equation 
is proved. See also [21], [11], [10] and [23]. However, since we cannot add 
constants to solutions, the Harnack estimates do not directly imply the local 
Holder continuity. To show that our proof is based on a general principle, we 
consider the case where the Lebesgue measure is replaced by a more general 
Borel measure, which is merely assumed to be doubling and to support a 
Poincare inequality. In the weighted case, parabolic equations have earlier 
been studied in [1], [2] and [20]. See also [8j. 

This kind of doubly nonlinear equations have been considered by Vespri 
[23| . Porzio and Vespri |19) . and Ivanov [14) . [15) . The known regularity 
proofs are based on the method of intrinsic scaling, originally introduced 
by DiBenedetto, and they seem to depend highly on the particular form 
of the equation. However, the passage from one equation to another is 
not completely clear. For other parabolic equations, the problem has been 
studied at length, see [3], [3], [7] and [22], and the references therein. 

The difficulty with equation (jl.ip is that there is a certain kind of di- 
chotomy in its behavior. Correspondingly, the proof has been divided in 
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1. Introduction 




2 < p < oo. 
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two complementary cases: 

Case I : < essinf n << essoscu 

and 

Case II : vF~'^ut ~ Cut. 

In large scales, i.e., in Case I, the scaling property of the equation domi- 
nates and, consequently, the reduction of the oscillation follows immediately 
from Harnack's inequality. In small scales, on the other hand, the oscillation 
is already very small and thus the solution itself is between two constants, 
the infimum and the supremum, whose difference is negligible. Correspond- 
ingly, the nonlinear time derivative term, which formally looks like vP~'^ut, 
behaves like ut and we end up with a p-parabolic type behavior. However, 
also in this case, we still need to modify the known arguments. In particular, 
the energy estimates are not available in the usual form and we need to use 
modified versions as in [5] , [T^ and [25] . 

Our argument also applies to doubly nonlinear equations of p-Laplacian 
type that are of the form 

-^^^ -V ■ A{x,t,u,Vu) = 0, 

with A{x,t, ■, •) satisfying the usual structure assumptions. For expository 
purposes, we only consider 

Very recently, a direct geometric method to obtain local Holder continuity 
for parabolic equations has been developed in [B] and [S]. Despite the effort, 
the general picture remains unclear. 



2. Preliminaries 

Let ;U be a Borel measure and Q be an open set in W^. The Sobolev space 
H^'P{^};fi) is defined to be the completion of C°°{^1) with respect to the 
Sobolev norm 

\\u\\i,P,n= (^j^{\u\P + \Vundf?j '\ 

A function u belongs to the local Sobolev space if it belongs 

to H^'P[Q'; fi) for every Q.' ^ Q. Moreover, the Sobolev space with zero 
boundary values Hq'^{Q,; fj.) is defined as the completion of C^{^1) with 
respect to the Sobolev norm. For more properties of Sobolev spaces, see e.g. 

m- 

Let tl < t2- The parabolic Sobolev space U'{ti,t2\H^'^{^; n)) is the 
space of functions u{x,t) such that, for almost every t, with ti < t < t2, the 
function u{-,t) belongs to H^''p{Q; fi) and 

' / {\u\P +\Vu\P)du < oo, 

where we denote di^ = dfj,dt. 

The definition of the space Lf^^{ti,t2; H^^^{Q; fj,)) is analogous. 



LOCAL HOLDER CONTINUITY 3 

Definition 2.1. A function u € L^^^{ti,t2', H^^^{Q; fi)) is a weak solution of 
equation (jl.ip in x (ti,t2) if it satisfies the integral equality 

^*'y" (^\Vu\P-^Vu-V(l)-uP-^^^ du = (2.2) 

for every (p G C^^iQ x (ti,t2))- 

Next, we recall a few definitions and results from analysis on metric mea- 
sure spaces. The measure /x is doubling if there is a universal constant 
Dq > 1 such that 

n{B{x,2r)) <Dofi{B{x,r)), 
for every B(x,2r) C 0,. Here 

B{x, r) = {y eW^ : \y - x\ < r} 

denotes the standard open ball in M°'. Let < r < R < oo. A simple 
iteration of the doubling condition implies that 

f^iB{x,R)) ^^(RV' ^ 



li{B{x, r)) \ / ' 

where = log2 Dq . A doubling measure in also satisfies the following 
annular decay property. There exist constants < a < 1 and c > 1 such 
that 

fi{B{x, r) \ B{x, (1 - 6)r)) < c6''fi{B{x, r)), (2.3) 

for all B{x, r) C and < (5 < 1. 

The measure is said to support a weak (l,p)-Poincare inequality if there 
exist constants Pq > and r > 1 such that 

/ \u-UB{x,r)\d^^<Por\f |Vn|Pd/i 

JB{x,r) \J B{x,Tr) J 

for every u G H^^{VL; ^u) and B{x,Tr) C il. Here, we denote 
UB{x,r) = i udn = ——J — / udfi. 

^ ' JB{x,r) f^{B{x,r}) JB(x,r) 

The word weak refers to the constant r, that may be strictly greater than 
one. In with a doubling measure, the weak (l,p)-Poincare inequality with 
some r > 1 implies the (l,p)-Poincare inequality with r = 1, see Theorem 
3.4 in |12) . Hence, we may assume that r = 1. 

On the other hand, the weak (l,p)-Poincare inequality and the doubling 
condition imply a weak (k, p)-Sobolev-Poincare inequality with 

d^-p ^ (2.4) 

where is as above. In other words, Poincare and doubling imply the 
Sobolev inequalities. More precisely, there are constants C > and r' > 1 
such that 

/ \u-UBix,r)rdfiy^^ <Cr(-f \Vu\Pdf}\ , (2.5) 

•JB{x,r) ' \JB{x,T'r) J 
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for every B{x,T'r) £ Q. The constant C depends only on p, Dq and Pq. 
For the proof, we refer to [12]. Again, by Theorem 3.4 in [12] we may take 



For Sobolev functions with the zero boundary values, we have the fol- 
lowing version of Sobolev's inequality. Suppose that u G H^'^{B{x,r)] ^). 
Then 

(\ i/« / \ i/p 

/ H^dA <Cr[i . (2.6) 

JB(x,r) J \JB{x,r) J 



For the proof we refer, for example, to 

Moreover, by a recent result in [TB], the weak (l,p)-Poincare inequality 
and the doubling condition also imply the (1, g)-Poincare inequality for some 
q < p, that is 

f \u- \dn<Cri-l I Vul'' dfj, ) . (2.7) 

JB{x,r) ' \JB{x,r) J 

Consequently, also (j2.5p holds with p replaced by q. We also obtain the 
{q, g)-Poincare inequality for some q < p. 

In the sequel, we shall refer to data as the set of a priori constants p, d, 
Dq, and Pq- 

Our main result is the following theorem. 

Theorem 2.8. Let 2 < p < oo and assume that the measure is doubling, 
supports a weak (l,p)-Poincare inequality and is non-trivial in the sense that 
the measure of every non-empty open set is strictly positive and the measure 
of every bounded set is finite. Moreover, let u > be a weak solution of 
equation (jl.ip in W^. Then u is locally Holder continuous. 

We will use the following notation for balls and cylinders, respectively: 

B{r) = B{0,r) 

and 

Qt{s,r) = B{r) x {t-s,t). 
For simplicity, we will also denote 

Q{s,r) = Qo{s,r) = B{r) x (-s,0). 
Recall Harnack's inequality from |[17j. 

Theorem 2.9. Let 1 < p < oo and suppose that the measure fi is dou- 
bling and supports a weak {l,p)-Poincare inequality. Moreover, let u > 
be a weak solution to (jl.ip in M"^. Then there exists a constant Hq = 
Ho{p,d,Do,PQ, {t-{s- rP))/rP) > 2 such that 

ess sup u < Hq essinf u, 

Qt{rP,r) Qs{rP,r) 

where s > t + . 

Proof See [n]. □ 

In addition, in [l?j it is also proved that all solutions of equation (jl.ip are 
locally bounded. In the sequel, we will assume this knowledge without any 
further comments. 
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3. Constructing the setting 

Our proof is based on the known classical argument of reducing the os- 
cillation, see [Ij, [7] and [22]. However, the equation under study has some 
intrinsic properties which are not present, for instance, in the case of the 
p-parabolic equation. In large scales, the scaling property dominates and 
the oscillation reduction follows easily from Harnack's inequality. In this 
case, the equation resembles the usual heat equation. 

In small scales, in turn, the equation changes its behavior to look more 
like the evolution p-Laplace equation. Indeed, when we zoom in by reducing 
the oscillation, the infimum and the supremum get closer and closer to each 
other. Consequently, the weight n^"^ in the time derivative term starts to 
behave like a constant coefficient and we end up with a p-parabolic type 
behavior. Resembling this divide between large and small scales, the proof 
has to be divided in two cases. 

We study the (local) Holder continuity in a compact set K and we choose 
the following numbers accordingly. Let 

/^o — '^-^ ^ /"d — ^^P ^' 

K K 

and define 

Furthermore, choose small enough so that 

{2Ho + 1)^0 < ^0 (3.1) 
holds. We will construct an increasing sequence {fJ-^} and a decreasing 
sequence {nf} such that 

fif - =uji = a'ujQ 
for some < cj < 1. Moreover, these sequences can be chosen so that 

ess sup u < fif 

and 

ess inf u > fi~ , 

for some suitable sequence {Q^} of cubes. Consequently, 

essoscn < Wj. 

The cubes here will be chosen so that their size decreases in a controllable 
way, from which we can deduce the Holder continuity. Observe also that if 

{2Ho + l)/iT < ujj, 
fails for some jo, the above sequences have been chosen so that 

^ <2Ho + 2 (3.2) 

for all j > jo. 

We are studying the local Holder continuity in a compact set K. Our aim 
is to show that the oscillation around any point in K reduces whenever we 
suitably decrease the size of the set where the oscillation is studied. 
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The next step is to iterate this reduction process. We end up with a 
sequence of cyhnders Q*. For ah purposes, in the sequel, it is enough to 
study the cyhnder := Q{rjrP, r) instead of the set K. Indeed, for any 
point in K we can build the sequence of suitable cylinders, but since we can 
always translate the equation, we can, without loss of generality, restrict the 
study to the origin. 

The equation (jl.ip has its own time geometry too, that we need to respect 
in the arguments. This is important, especially in small scales, when the 
equation resembles the evolution p-Laplace equation. We will use a scaling 
factor 7] = 2^i(P~^)+^ in the time direction, where Ai > 1 is an a priori 
constant to be determined later. 



We start the proof of Theorem 12.81 by proving the usual energy estimate 
in a slightly modified setting, which overcomes the problem that we cannot 
add constants to solutions, see [5], [15] and [25]. We introduce the auxiliary 
function 



4. Fundamental estimates 



J{{u-k)±) 





(4.1) 



<i-l{k + {u-kur-\u-k)i 

< ^-^u^~\u-k)l 



(4.2) 




(4.3) 



Observe that the assumption p > 2 is used here. 
On the other hand. 




(4.4) 



> 



(p-1) 



uP~\u-k)l. 
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Moreover, 

(u-k)- 



p(u—k)- 

j{{u-k).) = {p-i) / (A;-er2^d^ 

JO 



< (p- / (4.5) 





= (p-l)F — ^ — . 

Now we are ready for the fundamental energy estimate. 

Lemma 4.6. Let u > be a weak solution of (jl.ip and let k > 0. Then 
there exists a constant C = C{p) > such that 

ess sup / J{{u—k)±)(p^ dfi + / / \V{u—k)±ipfdi' 
ti<t<t2 Jn Jti Jn 

<c r [ {u-k)i\v^\Pdu+c r [ j{{u-k)±)^p-' (^) du, 

Jti Jn Jti Jn VCr/ + 



for every nonnegative ip € CQ"(r2 x {ti,t2)). 

Proof. Let ti < ri < T2 < t2- We formally substitute the test function 
(j) = ±{u—k)j-(p^ in the equation and obtain 

0=^"Y (^|Vnr^V^-V(/.+ ^^""J'^V ^ du 

= r I \V{u-k)\^-'^{±V{u-k)±)-V{±{u-k)±<fP)du (4.7) 

Jti Jn 

'■r2 r Ql^uP-^^ 



n Jn 



dt 



{u—k)±'f'P du. 



Now the first term on the right-hand side can be estimated pointwise from 
below as 

\V{u-k)\P~^{±V{u-k)±) ■ V{±{u-k)±^) 

> \V{u-k)±\PipP - p\V{u-k)±\P-^ipP-\u-k)±\Dip\, 

and the last term is estimated further by Young's inequality as 

- p\V{u-k)±\P~^ipP'\u-k)±\D^\ 

> -^|V(n-fc)±|V - C{u-k)l\D^\P. 



Thus we have 

^ r [ \V{u-k)±ip\Pdu± r [ ^^^^ ^\ u-k)±ipPdu 

Jti Jn Jti Jn 

<c r [ {u-k)i\v^\Pdu. 

Jt^ Jn 



2 



(4.8) 
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Using (j4.ip and integrating by parts, 



Ti Jn 9t J^^ Jq dt 

J{{u{x,t)-k)±)ipP{x,t) dn 

^ J t=Tl 

p r [ j{{u-k)±)^p-^^dv. 



Tl JO. 



So we obtain 



f J{iu{x,T2)-k)±)ipPix,T2)dfi+ r [ \V{u-k)±ip\Pdiy 
Jn Jti Jn 



<C J{{u{x,n)-k)±)ipP{x,n)dfi + C {u-k)l\V^\Pdiy (4.9) 

Jn Jti Jn 



+ C £ I^J{{u-k)±)^P-'(^^) du. 



Now we can drop the second term from the left hand side, let ri — > ti, choose 
T2 such that 

/ J{{u{x,T2)-k)±)^{x,T2)dii>]-esssn^ J{{u-k)±)^ d^i (4.10) 
Jn ^ ti<t<t2 Jn 

and estimate the limits of integration on the right hand side of (j4.9p . On 
the other hand, we can also drop the first term on the left hand side of (j4.9p 
and let ri ti and T2 — )• t2- Summing the estimates for both terms gives 
the claim. □ 

Let us denote 

few := *(fff c) = (in (_^^t_))^ . 

The following logarithmic lemma is used in forwarding information in time. 

Lemma 4.11. Let u > be a weak solution of equation (jl.ip . Then there 
exists a constant C = C{p) > such that 

esssup / uP-^'il;l{u)(pP dfi < [ kP-^iljl{u){x,ti)(fP{x) du 
ti<t<t2 Jn Jn 



+ 0^1 ij.\{ij.y\^-p\Vip\pdu 

Ju Jn 



and 



esssup / kP-'^^l{u)<fP dii< / uP-'^^l{u){x,ti)^P{x)dii 
ti<t<t2 Jn Jn 

+ 0^1 i>+\{i^+y\^-p\Vip\pdu. 

Jti Jn 

Above, (f G Cq°{Q) is any nonnegative time-independent test function. 



LOCAL HOLDER CONTINUITY 



Proof. Choose 

in the definition of weak solution and observe that 

(^|)" = (i + ^±)(V'i)2. 

The parabohc term will take the form 



(4.12) 



*2 /■ a_ 

dt 




ti Jn 



t2 f 

Ft 




kP- 



t2 



(f)±(r)r^ dr dfi 



1 *2 



n Jk 



Now an integration by parts gives 

(P±{r)rP-^ dr = ^ {iljl{r))' r^-^ dripP 



Jk 



In the plus case, we have 

(f)^{r)rP~'^ dr 

> 'tPl{u)uP~^ipP - V'i(n)(nP-2 - kP-^)ipP 
= {p-l)i^l{u)kP-^ipP 

and trivially 

^+{r)rP-^dr <4^l{u)uP-^ipP, 



since p > 2. Similar estimates are true also for the minus case. 

On the other hand, by using (|4.12|) together with Young's inequality, we 
obtain 

\Vu\P-^Vu ■ V<P± = \Vu\P-^Vu ■ V((V±(n))V^) 

> \Vu\P{l + ^±)(^±)V^ - 2p|Vn|P-V±V'±93P"V</'l 



1 



> -\vu\p{i + v±)(V'±)>^ - Ci,±{i:'^Y~p\Vip 

almost everywhere, from which the claim follows. 



n 



We will need the following notations in the next lemma, which is the most 
crucial part of the argument. Let 



Q± = S„ X = B{rn) X {t* - ^^rP,t 



and 



= {{x,t) e Qt : ±u{x,t) > ±kt} , 
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for n = 0, 1, 2, ... . 

Recall the definitions fif — jjL~ = coi = a^ujo, where fif > esssupg,: u and 
fj,~ < essinfgi u, for i > 1. Observe, however, that we have to choose 

IJ-o ^ ess sup It and //q > essinfit, 

K K 

where the infimum and the supremum are taken over K instead of . This 
is because we need the argument to be independent of the initial cylinder 
. Now we are ready to prove the fundamental lemma everything depends 
upon. 

Lemma 4.13. Let < ej- < 1, {k^)n be an increasing sequence and {k~)n 
a decreasing sequence, both of nonnegative real numbers. Suppose u > is 
a weak solution of equation (jl.ip . 

(u - kt)± < £±uJi and \k^+i - k^\ > 

In addition, assume further 

u > ^k- (4.14) 

and 

/i+<2A:+, n = l,2,... (4.15) 

for the minus and plus cases, respectively. Then there exist constants 
C_ = C{Do,Po,Co,p) > and C+ = C{Dq.Pq,p) > such that 

for every n = 0, 1, 2, . . . . Here k is the Sobolev exponent as in (j2.4p and 

Proof. Choose the cutoff functions (p^ € C^{Q^) so that < (p^ < 1, 
V'n = 1 in Qn+i and 



\V^t\<—^ and 



dt 



< ^ . (4.17) 



Denote in short 

"Vn — iu — kn)±, kji = kj^ , £ = £± 

and 

Qn = Bn xTn = Q^, An = , 7 = 7=^, Lfn = ft- 
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By Holder's inequality, together with the Sobolev inequality (|2.6p . we 
obtain 



J Qn + l 



dv 



_KQn)__/' 2{1-p/k)+p p{1-p/k)+p ] 



(4.18) 



<Ci (i vl^ldA U iVnifnTdl,) dt 

JT-n \J Bn / \J B„ J 

/ r \1-p/k /■ 

<CrP I esssupi vlipPdnj i \V{vnipn)\'' du. 
Here, we applied the doubling property of the measure v giving 

l^{Qn+l) ~ 

We continue by studying the term involving the essential supremum. By 
the assumption 

u > —K 

and ()4.4p . we obtain 

{u - K)l < -^u^-Pj{{u - k-)^) < C{k-f-Pj{{u - k-)^). 
p — 1 



On the other hand, the lower bound (|4.3p gives immediately 

{u-kt,)l<C{k+f-'^J{{u-kt)+). 
Using these estimates together with the energy estimate, Lemma 14.61 yields 

J Bn Tfi J Bn 

Furthermore, the estimates ()4.2p and ()4.5p imply 

J{{u-kn) + ) < C{ktr-Hu-kn)l, J{{u-kn)-) < C {k-y-\u-knf.. 

For the plus case, we used (|4.15p . Next, using (|4.17p . we arrive at 

-1 ( d(Pr, 



jQn V 



dt 



dv 



<C2-4 U+^^^vl\dv (4.19) 
1 / eoJi\^~^\ u{An) 



< C2''''{£uj,Y 1 + 



1 \ kn J / v{Qr 



where the last inequality follows from the fact that {u—kn)± < £±0Ji. Thus, 
we conclude 



ess sup 

Tn 



£.&f..,.,C2-(...f(.(|^)'-%i);4|;). (4.0) 
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Furthermore, since 

/ \y{VnVn)Tdv<cj \VVn\''^ldv + cj vl\V^n\Pdu, 

^ Qn Qn ^ Qn 

applying again the energy estimate and (|4.19p leads to 



/ 

JQ, 



|VK^„)rdz.<C2"P(e^,r ( 1 + i '] (4.21) 



To finish the proof, note first that 

>\h^ — 

>2-("+2)£^^.. 

It then follows that 



/ 

JQ 



,2(1-./.)+. > (2-i--^^)euy^'-'^^^~'' (4.22) 



n + l 



u{Q 



n+l 



Inserting estimates ()4.20p . ()4.2ip . and (|4.22p into (j4.18p concludes the proof. 

□ 

Remark 4.23. If we have the extra knowledge that {u — k~)- = 0, almost 
everywhere in B{r) at a given time level, we can choose the test functions to 
be independent of time and the cylinder Q~ so that the length in the time 
direction stays constant, and the bottom of the cylinder stays at the given 
time level. In this case, by choosing a time independent test function, the 
right hand side of the energy estimate simplifies so that we can get rid of 
the term +1 in the formulation and 

/ N p-2\ i-P/'^ 



in (|4.16p . This will get us the required extra room in the end of the first 

alternative of Case II. 

Furthermore, in the previous lemma, we chose the radii of the cylinder as 

r r 
+ 



2 2"+i 

However, the factor 2 in the denominator can naturally be replaced by any 
greater number. 

We start the proof by considering Case I. Here, we use the previous lemma 
only in the plus case. Consequently, the first case does not depend on the 
constant Cq. 

We recall a lemma on the fast geometric convergence of sequences from 

a- 

Lemma 4.24. Let (l^)n be a sequence of positive numbers, satisfying 

Yn+i < Cb^Y^+'^ (4.25) 

where C,b > 1 and a > 0. Then (l^)n converges to zero as n oo provided 

yo<C-W-"'. (4.26) 
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On several occasions in the sequel, we use this lemma, together with the 
fundamental estimate Lemma 14.131 to conclude that a ratio of the form 

converges to zero and consequently that v{A^) — )• as n — > oo. This will 
ultimately lead to a reduction of the oscillation which is our final goal. 

Once a recursive inequality of type ()4.16p has been established, the con- 
vergence to zero of i'{Af^) follows from the condition 

with 

a± = r;V(i-P/'^)^-i/(i-P/«)+i-{i-P/«)^^ (4.27) 

where the constants C± and T± are the constants from the previous lemma. 
Note that an explicit value of only follows after fixing C± and T±. 

5. The Case I 

Now we assume that ()3.ip holds. Our aim is to show that the measures 
of certain distribution sets tend to zero and that the local Holder continuity 
follows from this. 

We start by studying the subcylinder Q{r^,r) C Q{r]rf,r). Let 7^ = 1, 
£± = and 

, + _ + _ ^ _ ^0 

'^n — fJ-O 4 2n+2 ■ 

Observe that, after fixing these quantities, the constant can be fixed as 
well. 

We will study two different alternatives which are considered in the fol- 
lowing two lemmata, respectively. 

Lemma 5.1. Let A2 > 1 &e sufficiently large and letu >0 be a weak solution 
of equation (jl.ip . Furthermore, assume 

V ({(x, t) G B{r) X (-r^ -£) : u(x, t) > + = 0. (5.2) 

Then there exists a constant a E (0, 1) such that 

essosc u < aujQ. 

Proof. By the choices preceding the statement of this lemma, we have 

{u - k^)+ < e+LVQ. 
The assumption (jS.ip implies 

/^o = ^0 +^0 < 2t^o- 

Thus 

1 < < 4. 



14 KUUSI, SILJANDER AND URBANO 

Plug these in Lemma 14.131 to deduce 

On the other hand, by ()5.2p we have the trivial estimate 
choosing A2 > 1 sufficiently large. By Lemma 14.241 we conclude that 



as n — )■ 00. This implies 







ess sup u < fiQ — . 

Q((ir'i) ^ 

So, if this alternative occurs, we choose 

Wo 



and 



These choices yield 



as required, with 



^1 =1^0 



essosc u < [ 1 — - ] ojQ 



Qii 



2) '2/ 



3 

^=4- 



n 



For the second possibility, we have the following lemma. 
Lemma 5.3. Let u >0 be a weak solution of equation (jl.ip and suppose 

V ({(x, t) € B[t) X (-r^ -£) : u[x, t) > ^Uq" + > 0. (5.4) 

T/ien there exists a constant a = (t{Hq) G (0, 1) such that 

essosc u < (TWO- 



Proof. By assumption ()5.4p . we have 



. - Wo 
ess sup u > fiQ H . 

B(r)x(-rP,-^) 2 
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Now we can use Harnack's inequality (Theorem 12. 9p . together with the Case 
I assumption p.ip . to deduce 

essinf u > — ess sup u 



Hq 2Hq 
> /^Q H — 1 

Hq 2Hq + 1 2Hq 

too 



> 



2Ho{2Ho + l)- 

Observe that the constant Hq depends on A2, but this does not matter since 
A2 depends only on the data. 

Now, if we end up in this alternative, we choose 

- _ - wq 

^1 -^'^^2Hoi2Ho + l) 

and 

We also obtain 

Wo 

essosc U < , = crwn, 

" 2//o(2/^o + l) 

with 

1 

a = 1 



2Ho{2Ho + l)' 

as required. □ 



6. The case II 

In Case II the equation looks like the evolution p-Laplace equation. In 
this case, we need to use the scaling factor rj in the time geometry of our 
cylinders. The difficulty is now that we cannot use the Harnack principle 
anymore, as the lower bound it gives might be trivial. Indeed, the infimum 
can be larger than the lower bound Harnack's inequality gives. On the other 
hand, we have the following kind of elliptic Harnack's inequality. 

Suppose that jo is the first index for which assumption ()3.ip does not 
hold. Then we have 

</x+ < (2i?o + 2)/i-. (6.1) 

Clearly, this Harnack's inequality is valid also for every subset of the initial 
cylinder Qj^ = Q(rjrP,r) and, consequently, for every j > Jq. 
Recall, that = awjo-i and 



< < ^,„_i + (1 - (Tp, 



(2Ho + 2) - "^^0 - "^JO" , , 

2-a ^^-^^ 

< (2 - fT)Wjo_l < 



a 
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Thus, we obtain 



< ^—-^ < 1 



{2Ho + 2){2-a) " 2 - a coj^ 
and, consequently. 



- \ p-2 



where Ci = ^/(2 — aY ^- We will consider the cylinder 



By the above calculation, we have shrunk the cylinder by a factor which is 
controllable by the data. 

In the sequel, we will denote 

g = Cii^] . (6.3) 
\^jo J 

Recall the definitions 

Qt = Qt* {l^ri, r„) = BnXTn = B{rn) x {f - l^r^^t*) 

and 

At = {{x,t) e Qt : ±u > ±kt}. 

Now the proof will follow the classical argument of DiBenedetto, see [3] 
and ^22j, and is again divided into two alternatives. In the first one, we 
assume that there is a suitable cylinder for which the set where u is close 
to its infimum is very small. In the second alternative, we assume that this 
does not hold true. 

6.1. The First Alternative. We first suppose that there exists a constant 
oq € (0, 1) (to be determined in the course of the next lemma, depending 
only on the data) such that 

({(x,t) G Qo : n < + ^^}) < aou{Q^), 

for a cylinder 

Q- = Qt*{erP,r) C Q(2^i(P-2)+i0r^ r). 

Our aim is to use Lemma 14.131 to conclude for the reduction of the oscil- 
lation. 

Lemma 6.4. For every s > 3, 

.({(x,0GQ(^(9^9:n(x,t)<^-+^}) 



(s-3)2 V \4J '4 
where 9 is as in (16.31). 
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Proof. We start by using Lemma 14.131 with the choices 

£_ = 1/2 and 7~ = 6. We also need the assumption (|6.ip to deduce that 

u > ^(/xT^ + 2(2i/o + 2)fiJ^) > ^(/.T^ + 2u.,J > 1-k-, (6.5) 

with Co = 3{2Hq + 2). This knowledge is needed in Lemma [4. 131 Now, after 
fixing e_, 7~, A:~ and Co we can fix Oq , see ()4.27p . 
We also obtain, using ()6.ip and ()6.3p . the bounds 



2P-2 1 / 7^- \P-2 2P-2 

Ci 7 v^^-^ioy 

and thus we can conclude 

By the assumption of this alternative, together with the lemma of fast 
geometric convergence (Lemma I4.24p . we have u> k almost everywhere in 
Qt*i9{r/2)P,r/2). Thus 

(m-/c)_ = 

and consequently 

iPiu) := fin ( — ^ ^ =0 

almost everywhere in Qt* {9{r /2)p , r /2) . Let 

to < -e (-] (6.6) 

be a time level such that this is true for almost every x E B{r/2). 
Now our goal is to apply Lemma |4. 1 1 1 with 

k = ^'+^, c = ^ 

and 

= ess sup {u — A;)_, 
Q 

where Q = Q{r]9rP,r). Choose (p G Cg°(i?(r/2)) independent of time such 
that < (p < 1, if = 1 in B{r/A) and 

iVc^l < -. 
r 

In the set {u < /i~ + we have 

ip^>is- Sf W 2, 

and, on the other hand, 

V'<(s-2)ln2 and IV'?'^ < (^^''"^ 



18 KUUSI, SILJANDER AND URBANO 

The use of these estimates in Lemma 14.111 gives 

{f^-r~\s - 3)2 In^ 2 . f,{{x G Bir/A) : n{x, t) < + ^}) 

< ess sup / uP~'^^'^{u){x,t)Lp'^{x)d^ 

to<t<0 JB{r/2) 

< I kP-^^^{u){x,to)v>^ix)dfi 

JS(r/2) 

+ C f I il,\i;'\^-P\V^\Pd^idt 

Jto JB{r/2) 

<C(.-2)ln2(^^^J Mi?W4)) 

<C7(.-2)ln2( J^lj /i(i?(r/4)), 

for almost every t G (to) 0). Observe that, in the third inequahty, we plugged 
in r/ = 2''*i(^'~2)+^. The claim follows by integrating this estimate over 
(-0(r/4)P,O). 

□ 

We conclude this alternative with the following two lemmata. 

Lemma 6.7. Let u > be a weak solution of equation (jl.l|) and as- 
sume ()6.ip holds. Then 

u> u:; -\ TT a.e. in QiOl- 

where s depends only upon the data, and 9 is as in ()6.3p . 

Proof. Let 



8 2"+3' 
Q- = Bn X T = B{rn) x {to,0), 
where to is, as in the previous lemma, such that 

'r\P 



(u-/c)_(x,to) = 0, to<-e(£j^ 



for a.e. x G B{r/2). Moreover, define 



In this case, we obtain 

(u - k~)^ < e^Ujo, where e_ = ^. 

Observe also that 7" = -to/rP < r/6l = 2^i(p-2)+i^_ 

We will substitute these in Lemma 14.131 and, taking into account Re- 
mark S23] and estimate ()6.5p . we conclude as before that 
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Now choose s > Ai. Then, by Lemma [4.241 we have v{A~)/v{Q^) — )• as 
n — )■ oo, provided v{Aq)/i'[Qq) is smah enough. On the other hand, by 
choosing s large enough. Lemma [67il guarantees that i'{Aq)/v{Qq) can be 
chosen to be as small as we please. 
This gives 

and hence the lemma is proved. 

□ 

Lemma 6.8. There exists < cr < 1, depending only upon the data, such 
that 

essosc u < aoj^n- 



Proof. By the previous lemma. 



ess inf u > u- + , , , 

for some s > 1, which depends only upon the data and Ai. Observe that 
here we used the knowledge 



If this alternative occurs, we again choose 



and 

Finally, we get 
as required, with 



_ - ^jo 

1 

essosc n < I 1 



a = 1 



2^+1 
1 



2^+1- 



□ 



Remark 6.9. Here the choice of s is possible only after Ai has been deter- 
mined in the second alternative. Nevertheless, both of them are a priori 
constants which can be assigned explicit values depending only upon the 
data. 

6.2. The Second Alternative. In the second alternative, the assump- 
tion of the first alternative is not true. In this case, for every cylinder 
Qt*{6rP,r) C Q{ri6rP,r), we have 



u e Qt'{erP,r) : uix,t) > //T + ^} 



< (1 - ao). 



.10) 



u{Qt,{erP,r)) 

where ao := is the same constant as in the first alternative. This implies 
that, for every t* G ( — (?/ — l)0r^,O), there exists a time level to with 

t* - OrP <to<t* - —rP 
_ u _ 2 
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for which 

fi{{x G B{r) : u{x,to) > /cq"}) < ^fi{B{r)). -^^^ 

^ ~ Y 

Indeed, otherwise we would have 

{{{x,t) e Qt*{erP,r) : u{x,t) > k^}) 

> I /" ({a; G B{r) : u{x,t) > k^]) dt 

>{l-aoHQt'ierP,r)), 

which contradicts ()6.10p . 

This alternative is also based on Lemma 14.131 We choose Ai in the defi- 
nition of k^ large enough so that we can force i^{Aq) to be small compared 
to u{Q+). 

We start with forwarding the information of (j6.1ip in time. 

Lemma 6.12. There exists s* > 0, depending only upon the data, such that 

3ao 



^^ ({x G B{r) : u{x,t) > - ^}) < ^-—^^,{B{r)) 

2 



for almost all t G (to,0). 
Proof. Let 



— Mm 



2s+n ' Oo 2 

and 

= ess sup (u—k)^, 
Q 

where s and n will be chosen later and Q := Q{9r]rP,r). Our aim is again 
to use Lemma 14.111 to forward the information in time. We will need some 
estimates for doing this. 
Recall the definition 



V'+(u) = *(i^+,(u-A:)+,c)=ln^ 
Trivially, we have 



c + Ht-{u-ky 



V'+(n)<ln^ =nln2 



and, on the other hand, in the set 



we get 



The last estimate we need is 



- V2^ 
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Let now ip E CQ°{B{r)) be a cutoff function which is independent of time 
and has the properties 0<ip<l, {p = lm B({1 — 6)r) and 

|v*.| < J^, 

where < 5 < 1 is to be determined later. 

Apply Lemma 14.111 with these choices to conclude 

(n - if In^ 2 • fi{{x e B{{1 - 5)r) : u{x, t) > I}) 



<esssup / '4'^{u){x,t)ip'' {x) 

to<t<t* JBir) 
JBir) ^kJ 

+ Ck^-P f I ip+\{^+)'\^-P\V^\Pdfidt 

< ln2(2) (^^-A ^—^KB{r) x {to}) 
\ „+ ^ / 1 

for almost every t € {tQ,t*). Observe, that in the third inequality we 
used dSH]). 

Now, by the annular decay property (|2.3|) . we have 

fi{{x eB{r) : u{x,t) > I}) 

< n{B{r) \ B{{1 - 6)r)) + ^({x G B{{1 - 6)r) : u{x, t) > I}) 

< C5XB(r)) + fi{{x G 5((1 - 6)r) : n(x,t) > /}). 

For the first term, we choose 6 small enough so that 

ao 



< 



and for the second term we use the previous estimate. Indeed, by choosing 
s and n large enough so that 

1 — ao n / 1 \ 3ao 



< 1 



l_^{n-ir\l-±J - 8(l-^ 
and 

Cn f 1 Q-Q 



(ln2)5P(n- 1)2 ^2^ - 1 



< 



(-?)' 



16 



we get the claim for almost every t G {to^t*) with s* = s + n. Recall that 

t* - OrP <to<t* - —rP. 
_ u _ 2 
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Finally, since the above holds for every cylinder Qf {Or^, r) C Q{ri9rP, r), we 
can conclude that for almost every 

t > -rjOrP + 9rP - = (^-r/ + 1 - ^) Or^, 

we have 

1 _ 

G B{r) : u{x,t) > /}) < 4^f^{B{r)). 

~ T 

□ 

Remark 6.13. Now we can choose r] so large that the previous lemma holds 
for almost every t G ( — ^r^jO), i.e., 

— 71 +1 < 

and hence 

7? = 2^i(P-2)+i >2-ao. 
But this is always guaranteed for Ai > 1 and p >2. 

We are ready to prove the final estimate, which, together with Lemma r4.241 
gives the reduction of the oscillation. Let 

E,{t) = {xeB{r):uix,t)>f,l-'^} 



E, = {{x, t)eQ{ 4r', r) : u{x, t) > ^+ -^}. 



and 

Then we have the following lemma. 

Lemma 6.14. For every ai G (0, 1), there exists Ai > such that 



,y{Q{frP,r)) 



< ai. 



2 

Proof. Denote 
and 

k = 4-^, 

2^^ 

where s > will be chosen large. Let also 

h — k, u > h, 
V = ■{ u — k, k < u < h, 
0, u<k. 

By the previous lemma, we can choose s large enough, namely s > s*, so 
that, for almost every t G ( — ^r^,0), we have 

G -B(r) : t;(a;,t) = 0}) = ^({x G B{r) : ^(x,*) < A;}) 
> ^M-B(r)). 



LOCAL HOLDER CONTINUITY 
Thus, for almost every t £ 



23 



^r^, O) , we obtain 



VB{r) {t) = f vdfi< (l- ^ ) {h - k) 

JB{r)x{t} ^ 4 / 



and, consequently. 



"0 , 



h-k~ VB(r){t) > -^i^ - k)- 

Using the ((?, (7)-Poincare inequality for some q < p (see (|2.7p and the 
remark after that), yields 

{h-kyi^{Es+l{t)) < f— y / \v-VBir)itWdt, 
V«0/ JB{r)x{t} 

<Cr'i f \Vv\'idn = Cr'^ [ \Vu\'^ d^, 

JB{r)x{t} J Es(t)\Es+iit) 

for almost every t € ( — ^r^jO). The constant (4/ao)'^ above was absorbed 
into the constant C. Now we integrate the above inequality over time to get 

{h - kYv{Es+i) < Cr" I |Vu|« dv. 

J Es\Es+i 

Next, we introduce a cutoff function ip G C^{Q{7]6rP,2r)) such that < 
V'<l)V=liiiQ( ^rP,r] and 



C 

|V(/3| < — and 
r 



dip 



dt 



< 



Now Holder's inequality gives 
{h - kyv{Es+i) < Cr<i 



c 



q/p 



JEs\E,,j-i I 



g/p 



<CrM / \V{u-k)+\PipPdu] u{Es\E,+if~'ilP. 

IQ(r]erP,2r) J 



By choosing Ai > s > s* in the definition of r] large enough, the first 
factor on the right hand side can be estimated by Lemma 14.61 and (|4.2p as 



Ic 



\Viu-k)+\P^Pdi^ 



Q{rierP,2r) 

<C [ iu-k)l\Vip\Pdi^ 
+ C [ J{{u-k)+)ipP-' 

jQ{rierP,2r) 



dip 



dt 



dv 



.15) 



< — 



rP \ V 2 

C ^ujj,\p 
rP \ 2^* 



+ 



/ JQ(v9rP,2r) 



(n— /c)^ du 



24 KUUSI, SILJANDER AND URBANO 

In the last inequality we used the doubling property of the measure v. 
We obtain 

Finally, summing s over s*, . . . , Ai — 1 gives 

(Ai - s*)v{E^j'^^-'i^ < Cu {q (^rP^ryj"^" v {q (y^^r)) 



and hence 



Choosing Ai large enough finishes the proof. □ 

Lemma 6.16. Suppose that (|6.1U|) holds. Then there exists < o" < 1, 
depending only upon the data, such that 



essosc u < aoj 



Proof. Let 



JO- 



Q+ = i?(r„)x(-7+rP,0), 



r r 



2 2"+i 

and A+ as before. Substituting 7+ = 2^i(P-2)^^ ^ ]^/2Ai j^j^^j 



in Lemma [4. 131 and using (j6.ip to bound 

1< — <(2i7o + 2)P-2, 

7+ \e+ujjj 



yields 

< C 

By the previous Lemma, we can choose Ai large enough so that 

is as small as we please. Consequently, by Lemma [4.241 we obtain 

esssup n<^^.^-^, 

for some Ai > 1, which depends only upon the data. So if this alternative 
occurs, we choose 

+ + ^jo 

and 
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We also obtain 

essosc u < { 1 r — — - I Lv,„, 

as required, with 

a = 1 



2Ai+i- 

n 

Now the Holder continuity follows by standard iterative real analysis 
methods, see chapter III of [1], or [22]. 

Acknowledgments. The authors would like to thank Juha Kinnunen for 
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